No Maxwell Electromagnetic Wave Field Excited In Cloaked Concealment 
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The GL electromagnetic (EM) modeling is used to simulate the 3D EM full wave field propagation 
through cloaks. The 3D GL simulation of the EM wave field excited by a point source outside of the 
cloaks has been done. The simulation of the EM wave field from the point source excitation inside 
of the cloak device is presented in this paper. By using the GL modeling simulation, we found a 
phenomenon that there is no Maxwell EM wave field which is excited by nonzero local sources inside 
of the single layer cloaked concealment. The theoretical proof of the phenomenon by GL method is 
proposed in this paper. The GL method is fully different from the conventional methods. The GL 
method has double abilities of the theoretical analysis and numerical simulations to research the 
physical process and cloak metamaterial properties that is exhibited in this paper. 
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I. INTRODUCTION 



For the point source located outside of the cloak, the 
3D Global and Local field (GL) EM modeling method 
[1-3] has been used to simulate the 3D full EM wave 
field propagation through the cloaks. Our simulations 
and theoretical analysis in [1] verify the ideal cloak func- 
tions [4] . The EM wave field propagation outside of cloak 
does not penetrate into the concealment and never be 
disturbed by the cloaks. There are several papers to 
simulate the plane wave propagation through the cloak 
from outside of the cloak [5-7] . The plane wave is excited 
by plane source which can not be located inside of the 
cloak or concealment. The cloak simulation of the point 
sources and the local sources inside of the cloak is absent, 
in the mose published papers. The GL modeling simula- 
tions of the EM wave field through cloak and excited by 
the nonzero local sources inside of the cloak is presented 
in this paper. Moreover, by GL EM method simulation 
of the nonzero local sources inside of the concealment, 
we found a phenomenon that there exists no Maxwell 
EM wave field excited by nonzero local sources inside of 
the single layer cloaked concealment. The phenomenon 
is proved by the GL method theoretical analysis and the 
EM integral equation [1-3]. In paper [9] and [10], authors 
studied the effect on invisibility of active devices inside 
the cloaked region. Our statement is that "there exists 
no Maxwell EM wave field can be excited by nonzero 
local sources inside of the single layer cloaked conceal- 
ment with normal materials" . The detailed proof of the 
statement and 3D GL simulations are presented in this 
paper. 

Our GL method is fully different from conventional 
methods for cloak and physical and science simulations. 



It has advantages over the conventional methods. The 
GL method consistent combines the theoretical analyti- 
cal and numerical method together. In the GL modeling, 
there is no big matrix equation to solve and no absorp- 
tion condition on artificial boundary to truncate infinite 
domain. The method is a significant physical scattering 
process. The finite inhomogeneous domain is divided into 
a set of small sub domains. The interaction between the 
global field and anomalous material polarization field in 
the sub domain causes a local scattering wave field. The 
local scattering wave field updates the global wave field 
by an integral equation. Once all sub domains are scat- 
tered, the wave field in the inhomogeneous anomalous 
materials will be obtained. Therefore, the GL method 
can be used to both of theoretical analysis and numeri- 
cal simulation for physical and chemical phenomena and 
process. 

The description arrangement of the paper is as follows. 
The introduction is described in Section 1. The EM in- 
tegral equation is proposed in Section 2. The Global and 
Local EM field modeling is proposed in Section 3. The 
phenomenon that there is no the Maxwell EM wave field 
excited by the local sources inside of the concealment is 
proved in Section 4. The GL modeling simulations of the 
EM wave field propagation excited by the local source in- 
side of the cloak are presented in Section 5. We conclude 
this paper in Section 6. 
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II. 3D ELECTROMAGNETIC INTEGRAL 
EQUATION 



We have proposed the 3D EM integral equation in fre- 
quency domain in papers [1] and [2]. In this section, we 
propose the EM integral equation in time domain as fol- 
lows: 
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E(r, t) is the electric field, H (r, t) is the magnetic field, 
E b (r, t) and H b (r, t) is the incident electric and magnetic 
field in the background medium, E J (r' ,r,t) is electric 
Green's tensor, H J (r',r,t) is magnetic Green's tensor, 
they are excited by the point impulse current source, 
E M (r' ,r,t) and H M (r' ,r,t) are electric and magnetic 
Green's tensor, respectively, they are excited by the point 
impulse magnetic moment source, *t is convolution with 
respect to t, S [D] is the electromagnetic material param- 
eter variation matrix, 
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SDn and SD 22 are a 3 x 3 symmetry, inhomogeneous di- 
agonal matrix for the isotropic material, for anisotropic 
material, they are an inhomogeneous diagonal or full ma- 
trix, / is a 3 x 3 unit matrix, a(r) is the conductivity 
tensor, s(r) is the dielectric tensor, /2(r) is susceptibility 
tensor which can be dispersive parameters depend on the 
angular frequency lu, o b is the conductivity, e b is the per- 
mittivity, /Lth is the permeability in the background free 
space, n is the finite domain in which the parameter vari- 
ation matrix S [D] ^ 0, the (e(r) — e b I)E is the electric 
polarization, and (/2(r) — [i b I)H is the magnetization. 



III. 3D GL EM MODELING 

We propose the GL EM modeling based on the EM 
integral equation (1) in the time space domain. 

(3.1) The domain f2 is divided into a set of N sub 

JV 

domains, {fife}, such that O = (J ilk- The division can 

fe=i 

be mesh or meshless. 

(3.2) When k = 0, let E (r,t) and H (r,t) are the an- 
alytical global field, £ J (r', r, t), Hj(r', r, t), Eg* (r' , r, t), 
and Hq 1 (r 1 ,r,t) are the analytical global Green's ten- 



sor in the background medium, 
that Ek-i(r,t), H k -i(r,t), Ei 



By induction, suppose 



,t), 



E^L 1 (r f ,r,t), and H^ 1 (r',r,t) are calculated in the 



(fc — 1) step in the subdomain f2/c_i. 

(3.3) In {fife}, upon substituting 
JJ fc _i(r,t), E(_ l {r',r,t), 



Ek-i(r,t), 



E^ir'^t), 

and H^ 4 _ 1 (r' ,r,t) into the integral equation (1), the EM 



Green's tensor integral equation (1) in £1^, is reduced 
into 6x6 matrix equations. By solving the 6x6 matrix 
equations, we obtain the Green's tensor field E^(r', r, t), 
Hi (r' ,r,t), Ejf (r> ,r,t), and (r> ,r,t). 

(3.4) According to the integral equation (1), the elec- 
tromagnetic field Ek (r, t) and Hu (r, t) are updated by 
the interaction scattering field between the Green's ten- 
sor and local polarization and magnetization in the sub- 
domain Qk as follows, 
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(3.5) The steps (3.2) and (3.4) form a finite iteration, 
k = 1, 2, • • ■ ,N, the En (r, t) and Hn (t, t) are the elec- 
tromagnetic field of the GL modeling method. The GL 
electromagnetic field modeling in the time space domain 
is short named as GLT method. 

The GL EM modeling in the space frequency domain 
is proposed in the paper Q , we call the GL modeling in 
frequency domain as GLF method. 



IV. NO MAXWELL ELECTROMAGNETIC 
WAVE FIELD EXCITED IN CLOAKED 
CONCEALMENT 



Theorem: Suppose that a 3D anisotropic inhomoge- 
neous closed strip cloak domain separates the whole 3D 
space into three sub domains, one is the cloak domain 
f2 c with the cloak material; the second one is the cloaked 
concealment domain ftd with normal EM materials; other 
one is the free space outside of the cloak. If the Maxwell 
EM wave field excited by a point source outside of the 
concealment to be vanished in inside of the concealment, 
then there is no Maxwell EM wave field excited by the 
local sources inside of the cloaked concealment. 

The Maxwell EM wave field is the EM wave field which 
satisfies the Maxwell equation and continuous interface 
boundary conditions. We call the Maxwell EM wave field 
as the EM wave field and use inverse process to prove the 
theorem as follows: Suppose that there exists Maxwell 
EM wave field excited by the local sources inside the 
concealment with the normal materials, the wave field 
satisfies the Maxwell equation in the 3D whole space R 3 
which includes the anisotropic inhomogeneous cloak do- 
main Q c and concealment fid, and satisfies the contin- 
uous interface conditions on the inner boundary surface 
Si and outer boundary surface S2 of the cloak domain 

Let R c = R 3 - n c U fl d , R d = R 3 - ft d , and by the EM 
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integral equation (1), the EM wave field satisfies 
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where G^j^ (r', r, i) is the EM Green's tensor, its compo- 
nents E J ,H J , E M , and H M (r', r, t) arc the EM Green's 
function on fl c 1J fid 1J R c , excited by the point impulse 
sources outside of the concealment, r G Rd- By the 
assumptions, G^^(r',r, t) exists on fl c \J fid \J R c and 

G'e^h ( r '> r ^) = 0' when r' G fid- The integral equation 
(5) becomes to 
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We consider the Maxwell equation in Rd, the virtual 
source located r, r € Rd and the point source located r s , 
r s G fid and r s ^ Rd, we have 



Vx 



-Vx 



[D] G J d M H (r',r,t)+IS(r',r)S(t), 



G J d M H (r',r,t) 



(7) 



and 



-Vx 
= [D b ] 



Vx 

E b 
H b 



E b 



(r',r s ,t) 



(r',r s ,t) 



(8) 



By using [E b (r,t) , H b (r,t)] to convolute (7), and 
Ggj(r',r,f) to convolute (8), to subtract the second 
result equation from the first result equation and make 
their integral in fl c 1J R c , and make integral by part and 
some manipulations, we can prove 
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dS. 



Because G^ H (r',r,t) = 0, r 1 G Qd, by continuous 
interface conditions of G^^(r',r, t), the term in right 
hand side of (9) is vanished, we have 
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Upon substituting integral equation (10) into the inte- 
gral equation (6), we have 



E(r,t) 
H(r,t) 



0. 



(11) 



From the continuous property of the EM wave field, 
we obtain the following over vanish boundary condition 
on the boundary S\ of the concealment Vld, we have 



E(r,t) 
H(r,t) 



= 0. 



(12) 



Si 



Because the EM wave field is excited by local sources 
inside of the concealment domain f^, it satisfies the fol- 
lowing Maxwell equation, 
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(13) 



The equation (13) and its over vanish boundary 
condition (12) form contradiction equations on the 
concealment fid, where [D] = diag [e fi] {d / dt) with 
the normal EM material parameters e and fi, r s G fid 
is the source location, Q(r',r s ,t) is the nozero lo- 
cal source inside of fid- There is no any Maxwell 
EM wave field to satisfy the contradiction equations 
(13) in fid with over vanish boundary condition 
(12) on the boundary Si. Therefore, we proved that 
there is no Maxwell EM wave field excited by the 
nonzero local sources inside of the cloaked 
concealment. 

In paper [1], we have used the simulations and the- 
oretical analysis by the GL method to prove that the 
excited EM wave field by the point source outside of 
the sphere cloak [4] through the cloak and never propa- 
gate enter the concealment and never be disturbed by the 
cloak. The assumptions and conditions in the theorem in 
this paper are validated, therefore, we stress that there is 
no Maxwell EM wave field excited by the nonzero local 
sources inside of the center sphere concealment cloaked 
by sphere cloak and the arbitrary geometry closed strip 
cloak. In practical cloak metamaterial fabrication and 
experiments, the phenomenon is presented in this paper 
should be received attentions, because the EM wave field 
excited in the concealment may be an irregular EM chaos 
propagation in this region, it may interfere the EM de- 
vices and equipments working inside of the concealment, 
the high frequency irregular EM chaos radiation may hurt 
the health of the human. 



V. SIMULATION OF THE EM WAVE FIELD 
THROUGH THE CLOAK BY GL METHOD 

The simulation model: the 3D domain is 
-0.5m, 0.5to] x [-0.5m, 0.5m] x [—0.5™, 0.5m], the 
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Electric Wave Exx Propagation in Sphere EM Cloak 
In XY Plane Excited by Point Source Outside of the Cloak 
Using 3D GL Modeling in GLGEO 




Horizontal X (Meter), Made by GLGEO 

FIG. 1: (color online) By the point source excitation in x — 
1.1m on X axis outside of the cloak, at the time step llOdt, 
the electric wave field Exx propagation does disperse and split 
into the two phases around the concealment with an antenna 
0, and does not penetrate into it. 

Electric Wave Exx Propagation in Sphere EM Cloak Excited by 
Point Source in Cloak Using 3D GL Modeling in GLGEO 




Horizontal X (Meter), Made In GLGEO 

FIG. 2: (color online) The electric wave field Exx is excited 
by the point source inside of the cloak at x=0.25m on X axis, 
at the time step 15dt, the Exx propagates around the con- 
cealment and does not penetrate into it. 

mesh number is 201 x 201 x 201, the mesh size is 0.005m. 
The electric current point source is defined as 

S{r - r s )6(t)e, (14) 

where the r s denotes the location of the point source, 
the unit vector e is the polarization direction, the time 
step dt = 0.3333 x 10~ 10 second, the largest frequency 
/ = 10GHz, the shortest wave length is 0.03m. The EM 



cloak fi c is the spherical annular with the center in the 
origin and internal radius R\ — 0.2m and exterior radius 
R 2 = 0.3m. The cloak is divided into 50 x 50 x 50 cells, 
the antenna subdomain 0, f2 ff , is divided into 63 cells. 
The spherical coordinate is used in the sphere r < R2, the 
Cartesian rectangular coordinate is used in other where 
to mesh the domain. 

For a point source located outside of the cloak, the 
3D GL EM modeling has been used to simulate the EM 
wave field through the sphere, ellipsoid, cylinder, and 
arbitrary closed strip complex geometry cloaks, the sim- 
ulations and theorems of the single and multiple sphere 
cloaks are proposed in paper [1]. The cloak simulations 
in papers [5-7] are proposed for the outside plane wave 
through cloak. The plane source to excite the plane wave 
can not be located inside of the cloak and concealment. 
Simulation for the EM wave field excited by the point 
source, in particular, the source located inside of cloak, 
r s G il c , or inside of concealment r s 6 is lack. The 
Figure 1 shows that the electric wave field Exx is ex- 
cited by the current point source in the direction e = x 
and located at (1.1m, 0.0, 0.0) on the X axis where is 
outside of the cloak. At the time step llOdt, the elec- 
tric intensity wave field Exx is propagating through the 
sphere annular cloak and around the sphere concealment, 
it does disperse and split into the two phases around the 
sphere concealment, the front phase speed exceeds the 
light speed; the back phase is slower than the light speed. 
The wave front outside of the cloak is the same as the 
exact Exx propagation in free space, the Exx wave field 
outside of the cloak never been disturbed by the cloak 
and never penetrate enter the centre sphere concealment 
with the antenna ©. The simulations of the EM wave 
field through the cloak excited by the point source inside 
of the cloak are presented in this paper. The electric wave 
field Exx excited by the current point source in the di- 
rection e = x and located at the point (0.25m, 0,0) inside 
of the cloak is propagating at the time step 15<ii that is 
shown in Figure 2. The Figure 3 shows that the electric 
wave field Exx propagates around the cloaked conceal- 
ment at the time step 18dt , but does not penetrate into it 
and its right wave front has been outside of the cloak. At 
the time step 48dt, the electric wave field Exx is prop- 
agating outside of the cloak that is shown in Figure 4, 
however, to compare an Exx excited by the same source 
in free space, the electric wave field Exx propagation out- 
side of the cloak is disturbed by the cloak. Figures 2 - 
4 show that the wave field Exx propagates through the 
cloak and never enter the concealment. The concealment 
is complete concealed by the cloak from the EM wave 
field excited by the point sources inside of the cloak and 
in free space outside of the cloak. We did use the GL 
modeling to simulate many cases of the EM wave field 
excited by the point sources inside of the concealment. 
However, all simulations are unstable and chaos. When 
the EM wave field is propagating to arrive the interface 
boundary S± between the concealment and cloak, the GL 
simulation become unstable and chaos. The GL modeling 
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simulation experiments reminder us to think may there 
is no any Maxwell EM wave field excited by the point 
source inside of the concealment. This is the motivation 
and idea of our theorem proposed in this paper. After 
the rigorous proof by the GL method, we obtained the 
theorem that there is no Maxwell EM wave field excited 
by the point sources located inside of the concealment. 



Electric Wave Exx Propagation in Sphere EM Cloak Excited by 
Point Source in Cloak Using 3D GL Modeling in GLGEO 




Horizontal X (Meter), Made In GLGEO 

FIG. 3: (color online) At the time step 18di, the Exx propa- 
gates around the concealment and does not penetrate into it, 
its right wave front is outside of the cloak 



Electric Wave Exx Propagation in Sphere EM Cloak Excited by 
Point Source in Cloak Using 3D GL Modeling in GLGEO 
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FIG. 4: (color online) At the time step 48dt, the Exx is 
propagating outside of the cloak, the wave field propagation 
in the free space outside of the cloak is disturbed by the cloak. 

VI. CONCLUSIONS 

The GL method is used to simulate the invisibility of 
the sphere and arbitrary cloaks and theoretically and rig- 
orously proved theorem that there is no Maxwell EM 
wave field excited by the nonzero local sources inside of 
the concealment which is cloaked by the sphere cloak and 
arbitrary closed strip cloak. A least square or regulariz- 
ing chaos propagation of the EM wave field excited in 
the concealment will be modeling and inversion by the 
GL metre carlo method [8] in next paper. However, any 
field excited in the cloaked concealment can not be prop- 
agation outside of the concealment. 

The GL EM modeling is fully different from FEM 
and FD and Born Approximation methods and overcome 
their difficulties. There is no big matrix equation to 
solve in GL method. Moreover, it does not need arti- 
ficial boundary and absorption condition to truncate the 
infinite domain. The GL EM method consistent com- 
bines the analytical and numerical approaches together. 
The GL method has double abilities of the theoretical 
analysis and numerical simulations that is shown in this 
paper. 

The 3D GL simulations of the EM wave field through 
the single and multiple sphere, cylinder, ellipsoid, and 
arbitrary geometry cloaks show that the GLT and GLF 
EM modeling are accurate, stable and fast, ft saves more 
storages than the conventional methods and needs 10 to 
50 minute to run the 3D EM wave field through the cloaks 
with 64 to 128 frequencies in the PC. The high perfor- 
mance GL parallel algorithm in PC cluster and super 
parallel computer is very fast and powerful to simulate 
complex and large scale physical and chemical process. 

The 3D and 2D GL parallel software is made and 
patented by GLGEO. The GL modeling can be extended 
to its inversion [8[ and GL EM quantum field modeling 
to solve quantization scattering problem of the electro- 
magnetic field in the dispersive and loss metamaterials, 
cloaks and more wide anisotropic materials. 
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